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Abstract 

From e + e~ collision data acquired with the CLEO-c detector at CESR, we search for the non- 
DD decays ^(3770) — > jXcj, with Xcj reconstructed in four exclusive decays modes containing 
charged pions and kaons. We report the first observation of such decays for J = with a branching 
ratio of (0.73 ± 0.07 ± 0.06)%. The rates for different J are consistent with the expectations 
assuming ■0(3770) is predominantly a l 3 Di state of charmonium, but only if relativistic corrections 
are applied. 
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Observation of the narrow X(3872) and V(4260) states [l| above open charm threshold, 
and their possible interpretation as states beyond the traditional cc model of charmonium 
0, calls for thorough investigation of the lightest charmonium state above the DD threshold 
- -0(3770). The common interpretation of the -0(3770) assumes it is predominantly the l 3 Di 
cc state, with a small admixture of 2 3 S' 1 . Except for the large DD decay width and rough 
agreement with the potential model mass predictions, there have been no other experimental 
data to verify this assumption. Although decays of ■0(3770) to 7r + 7r~ J/-0, n°n 'J /-0 and r/J/tp 
have been measured to be non-zero 0, H| , such hadronic modes present a less sensitive probe 
of the charmonium model than rates for 0(3770) — > 7Xcj since they involve hadronization 
probabilities. 

Previously, we have reported observation of 0(3770) — > -fXa with Xci — * iJ/^k, J 'fat —> 
l + l~ The branching ratio for -0(3770) — > •yXco is predicted to be the largest p, IE El ; 
but the small branching ratio for Xco 7^/"0 reduces the sensitivity so much that only a 
loose upper limit could be set in Ref. p. However, hadronic Xco decays are copious and 
thereby offer complementary probes for these photon transitions. Backgrounds from DD 
decays and continuum processes are suppressed by full reconstruction of Xcj decays to a few 
exclusive hadronic final states. We use the following decay modes: Xcj K + K~ (2K), 
XcJ - * 7T + 7r~7r + 7r~ (47r), Xcj —* K + K~tt + tt~ (2K2n) and Xcj — > 7t + 7t~7t + 7t~7t + tt~ (6tt). To 
minimize sensitivity to large uncertainties in branching fractions and resonant substructure 
for these channels, we measure the rates relative to those seen in ip(2S) decays with the 
same detector, 

R ^(^(3770) -> 7Xcj) x B(xcj -j> ^K*) 
J B(^(2S) - 7 Xcj) x B( Xc j - tt±, K±) ' 

and normalize to B(ip(2S) — > 7Xcj) which was measured by fitting inclusive photon 
energy spectra. Thus, our results for $(-0(3770) — ► -fXcj) are not only independent of 
B(xcj 7t ~ ± ,K ± ), but also depend only on ratios of detection efficiencies for 0(3770) and 
ip{2S). The latter are almost independent of the resonant substructure and, therefore, can 
be more reliably determined. 

The data were acquired at a center-of-mass energy of 3773 MeV with the CLEO-c detector 
[ill ] operating at the Cornell Electron Storage Ring (CESR), and correspond to an integrated 
luminosity (number of resonant decays) of 281 pb _1 ((1.80 ±0.05) x 10 6 ) at the 0(3770) and 
2.9 pb" 1 ((1.51 ± 0.05) x 10 6 ) at the 0(25). The CLEO-c detector features a solid angle 
coverage of 93% for charged and neutral particles. The cesium iodide (Csl) calorimeter 
attains photon energy resolutions of 2.2% aX E 1 = 1 GeV and 5% at 100 MeV. For the 
data presented here, the charged particle tracking system operates in a 1.0 T magnetic field 
along the beam axis and achieves a momentum resolution of 0.6% at p = 1 GeV. Particle 
identification is performed using Ring-Imaging Cherenkov Detector (RICH) in combination 
with specific ionization loss (dE/dx) in the gaseous tracking volume. 

We select events with exactly 6, 4 or 2 charged tracks and at least one photon candidate 
with energy above 60 MeV. The highest energy photon is considered to be the signal photon, 
while other neutral clusters in the calorimeter are considered fragments of hadronic showers, 
and therefore ignored. We separate pions and kaons using a log-likelihood difference, which 
optimally combines the dE/dX and RICH information. The track is considered a kaon if the 
kaon hypothesis is more likely. The RICH information is used only if the track momentum 
is above kaon radiation threshold (700 MeV) and the number of Cherenkov photons for the 
kaon hypothesis is required to be at least 3 in this case. We also impose 3cr consistency 
on dE/dx. Those tracks not identified as kaons become pion candidates if they satisfy 3a 
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TABLE I: Efficiencies for ij){2S)/^} (3770) — > 7Xcj 5 Xcj — ► 7r=t )-^ =l= ) based on Monte Carlo of phase- 
space Xcj decays (i.e. no intermediate resonances). 





Efficiency (%) 
J=2J=1J=0 




IXcJ ~ 


-> 47T 


33 


35 


34 






2K2ir 


25 


27 


28 






6ir 


23 


25 


27 






2K 


43 


44 


42 


0(3770) - 


-> IXcJ - 


-> An 


35 


36 


34 






2K2it 


29 


30 


29 






6ir 


27 


28 


27 






2K 


44 


44 


41 



consistency with dE/dX. Events with odd numbers of kaons or pions are rejected. The total 
energy and Cartesian components of momentum of the selected charged particles and the 
photon must be consistent within ±30 MeV with the expected center-of-mass four-vector 
components, which take into account a small beam crossing angle. To improve resolution on 
the photon energy, we then constrain these quantities to the expected values via kinematic 
fitting of events. Selection efficiencies obtained with GEANT l3| based simulation of 
detector response are given in Table UJ 

The energy of the photon candidates is plotted for the data for different decay channels in 
Fig.n &n d Fig. |21 Fits used to extract signal amplitudes are also shown. Each photon line is 
represented by a detector response function, parameterized by the so-called Crystal Ball line 
(CBL) shape. CBL is a Gaussian (described by the peak energy, E , and energy resolution, 
ge) turning into a power law tail, 1/ (E —E+ const)™, at an energy of E —aaE- We fix a and 
n to the values determined from the signal Monte Carlo. The peak amplitude (A^^s)^)} 
peak energy and widths are free parameters in the fit to the ip(2S) data. The smooth 
background is represented by a first order polynomial. In the fit to the 0(3770) data only 
the peak amplitudes (A/>(3770)) are free parameters, while the CBL parameters are fixed to the 
predictions from the signal Monte Carlo. In addition to the smooth backgrounds, represented 
by a second order polynomial, the 0(3770) data also contain radiatively produced ip(2S) 
background. After our selection cuts, the latter cannot be distinguished from the ■0(3770) 
signal. They are explicitly represented in the fit by peaks with the amplitudes, A™^^ 71Q \ 

fixed to the values estimated from the ip(2S) data (A^^s)^) an d extrapolated to the -0(3770) 
beam energy with help of the theoretical formulae: 

A™(2S) = £</>(3770) • e</>(3770) • &x ' T ee {^{2S)) ■ I(s) 

Us) = / W(s, x) ■ b(s'(x)) ■ F x (s'(x))dx. 
Jo 

Here, we are using the same notation as in Ref. £ is the integrated luminosity; e is 
the efficiency; Bx is the branching ratio for ip(2S) ^Xcj — * lX {X is the hadronic final 
state) at the ip(2S) resonance peak; x is energy radiated in e + e~ — > ■yip(2S) divided by its 
maximal possible value (i.e. by -Ebeam = V^/2); s' is the mass-squared with which the ip(2S) 
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FIG. 1: Distribution of photon energy for 4.7T (top) and 2K2tt (bottom) decay samples in CLEO-c 
ip(2S) (left) and ^(3770) (right) data. Solid histogram is data, smooth curve is fit to the data. 
Dashed line shows radiative return background contribution from ip(2S) tail and dotted line is 
polynomial background. 



is produced (s'(x) = s(l — x)); W(s,x) is the initial state radiation probability (see Ref. 
for the definition and discussion); b(s') is the relativistic Breit-Wigner formula describing 
the ip{2S) resonance (b(s') = 12ttT r /[(s' - Ml) 2 + M|r|]); and F x (s') is the phase-space 
factor between the ip(2S) produced with a/s 7 mass and with its nominal mass, M R . F x (s') is 



equal [13[ to (i? 7 (s')/£' 7 (M|,)) 3 , where E 1 is the photon energy in il)(2S) — » 7Xcj decay. The 
?p(2S) nominal mass (Mr) and total width (r^) are taken from PDG 
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FIG. 2: Distribution of photon energy for 6tt (top) and 2K (bottom) decay samples in CLEO-c 
ip{2S) (left) and ^(3770) (right) data. Solid histogram is data, smooth curve is fit to the data. 
Dashed line shows radiative return background contribution from ip(2S) tail and dotted line is 
polynomial background. 

is taken from the CLEO determination utilizing e + e~ — > _yip(2S) at Eqm — 3773 MeV 
with ip(2S) decaying to J/ip through a hadronic transition [4|]. The radiative flux, W(s,x), 
strongly peaks for x — > making the ip{2S) background indistinguishable from the ?/>(3770) 
signal within our photon energy resolution. Unlike in our X = -yJ/ip analysis [Hf, where we 
used the published CLEO results for Bx and relied on the absolute value of the detection 
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TABLE II: Fitted signal yields for V(25)/i/>(3770) -> 7X c j, XcJ -> vr ± , J fC ± . The total number of 
the estimated i{)(2S) background events in the ^(3770) data {A™^^ 77 ^ ) is also given. The errors 
on the latter quantities are systematic. All other errors are statistical. 



Decay 

Mode J = 2 



Events 

J = 1 J = 



A 



in V(25) 
V(2S) 



.4 



in ^(3770) 
V(2S) 



^(3770) 



4?r 534 ±27 291 i 19 981 ± 36 
2K2ir 261 ±16 187 ± 14 745 ± 29 
6vr 469 ± 23 408 ± 21 744 ± 30 
2K 64 ±8 - 346 ±19 
All 1329 ± 40 886 ± 32 2816 ± 58 



All 25 ±6 12 ±3 25 ±6 



4tt 


9 ±10 


14 i 9 


112 i 16 


2K2lT 


6±8 


25 i 9 


73 i 14 


6tt 


5 ±12 


16 i 11 


65 i 16 


2K 


Oil 




24 i 6 


All 


20 i 18 


54il7 


274 i 27 



efficiency (e^(3770)), in this analysis we set 

4 in V>(2S) 

B x ~ V,(2S) 



ev(2S) ' ^V(25) ' 

where A^^s)^ is the signal yield in the fit to the ip(2S) data. Therefore, our estimates of 
the ip(2S) radiative tail background, 

.in V(3770) _ A in f(2S) e ^(3770) £y>(3770) p , ,,x 

^(25) iV^(2S) 

do not rely on absolute values of efficiencies, but only on their ratio between the -0(3770) and 
ip(2S) data samples. The upper range of integration in the definition of J(s) is a; cut ~ 30 
MeV/1887 MeV=0.016, because of our cuts on total energy and momentum. The signal 
yields in the ip(2S) and -^(3770) data are given in Table HTl 

The results for the ratio of branching ratios, Rj, for individual decay modes are given 
in Table IIHI Average values are calculated using inverse-of-statistical-errors-squared for 
weights. To estimate the statistical significance of -0(3770) — > ^Xcj signals, we fit the 
-0(3770) data with the background contribution alone and compare the fit likelihoods to our 
nominal fits. Combining likelihoods for all the channels, we obtain statistical significance 
of 1.3, 3.6 and 12.6 standard deviations for J = 2, 1 and 0, respectively. The sum of the 
photon spectra over the individual channels is shown for ip(2S) and -0(3770) data in Fig. 01 
Since no significant signal is observed for J = 2, we set an upper limit for this state. 

Various contributions to the systematic errors are listed in Table |TV| We simulated 
signal events assuming various resonant substructures and compared the efficiency ratio to 
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FIG. 3: Distribution of photon energy in CLEO-c ip(2S) (top) and -0(3770) (bottom) data summed 
over all analyzed modes (data points). The smooth curve shows the sum of the fits performed to 
the individual modes. The dashed curve shows the radiative tail from ip{2S). The dotted line 
shows the polynomial background. 



our nominal values obtained with the phase-space model to evaluate the error in efficiency 
simulation. Including the systematic errors, our results for the ratio of branching ratios 
are: R = (7.9 ± 0.8 ± 0.6)%, R x = (4.3 ± 1.6 ± 0.6)% and R 2 < 2.2% (90% C.L.). The 3% 
uncertainty in the number of ip{2S) resonant decays contributes to the Rj measurement, but 
cancels when multiplied by the inclusively measured B(ijj(2S) —> 7Xcj) (ifj- The results for 
^(^(3770) -> 7x c j) are (0.73±0.07±0.06)%, (0.39±0.14±0.06)% and < 0.20% (90% CL.) 
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TABLE III: The ratio Rj = £(^(3770) - 7Xc j, Xcj - vr±, K^/BtypS) - 7Xc j, Xcj - vr±, tf±). 
Only statistical errors are given here. 



Decay 


Rj in % 


mode 


J = 2 J = 1 J = 


47T 


1.3 ±1.5 3.8 ±2.6 9.6 ±1.4 


2K2tt 


1.7 ±2.4 9.9 ±4.0 8.2 ±1.7 


6tt 


0.7 ± 1.8 2.9 ±2.2 7.4 ± 1.8 


2K 


0.0 ± 1.4 - 6.0 ± 1.6 


Average 


0.8 ±0.8 4.3 ± 1.6 7.9 ±0.8 



TABLE IV: Systematic errors and their sources. 





Relative chanj 


j;e in % 




J = 


2 J= 1 


J = 


Luminosity 


1 


1 


1 


V>(3770) cross-section 


3 


3 


3 


Number of ip(2S) decays 


3 


3 


3 


Resonant substructure 


2 


< 1 


< 1 


±25% change in ip(2S) bkg. 


39 


6 


2 


Fit systematics 








±7% change in oe 


10 


8 


4 


±10% change in fit range 


17 


5 


1 


Using Gaussian signal shape 


9 


2 


1 


Decreasing bin-size to half 


15 


3 


< 1 


±1 order of bkg. polynomial 


47 


9 


2 


Total fit systematics 


53 


12 


5 


Total systematic error on Rj 


66 


14 


7 


B&(2S) -> 7Xcj) 


6 


5 


4 


Number of ip(2S) decays 


-3 


-3 


-3 


Total systematic error on 








£(^(3770) 7X C J) 


66 


15 


8 



for J = 0, 1 and 2, respectively. They are consistent with the results obtained previously 
by CLEO @ using X cj -» 7<//V decays: < 4.4% (90% C.L.), (0.28 ± 0.05 ± 0.04)% and 
< 0.09% (90% C.L.), correspondingly. The two analyses are complementary. While this 
analysis offers much better sensitivity for J = 0, the previous analysis is more sensitive for 
J = 1 and 2. The J = 1 signal is observed in both analyses. Combining both analyses we 
obtain 23(^(3770) -> jxa) = (0.29 ± 0.05 ± 0.04)%. 

We turn the branching ratio results to transition widths using r to t = (23.6 ± 2.7) MeV 
from PDG 14]. The results are given in Table IV| where they are compared to theoretical 



9 



TABLE V: Our measurements of the photon transitions widths (statistical and systematic errors) 
compared to theoretical predictions. The J = measurement comes from this analysis. The J = 2 
upper limit comes from Ref.jjJ. The J = 1 measurement comes from the combination of this 
analysis and of the result in Ref. [jj • 





r(V(3770) -> 7Xcj) m keV 
J = 2 J=l J = 


Out rpmilf^ 


< 21 


70 ± 17 


172 ± 30 


Rosner (non-relativistic) [7] 


24 ±4 


73 ±9 


523 ± 12 


Ding-Qin-Chao [6] 








non-relativistic 


3.6 


95 


312 


relativistic 


3.0 


72 


199 


Eichten-Lane-Quigg [8] 








non-relativistic 


3.2 


183 


254 


with coupled-channels corrections 


3.9 


59 


225 


Barnes-Godfrey-Swanson [9] 








non-relativistic 


4.9 


125 


403 


relativistic 


3.3 


77 


213 



predictions. 

The theoretical predictions are based on potential model calculations 0| of the electric 
dipole matrix element < l 3 Pj\r\l 3 Di >: 

where tQ is the c quark charge and a is the fine structure constant. The spin factors Cj 
are equal to 2/9, 1/6 and 1/90 for J = 0, 1 and 2, respectively ,15]. The phase-space 
factor (E 3 ) also favors the J = transition. Together, the spin and phase-space factors 
predict enhancement of the J = width by a factor of ~ 3.2 and ~ 85 over J = 1 and 
J = 2, respectively. In the non-relativistic limit, the matrix element is independent of J. 
The measured ratios of the widths, T /Ti = 2.5 ± 0.6 and r /r 2 > 8 (90% C.L.), are 
consistent with these crude predictions, therefore, providing further evidence that ?/>(3770) 
is predominantly a 1 3 D\ state. A small admixture of 2 3 S\ wave, necessary to explain the 
observed T ee (■0(3770)), is expected to increase T and T 2 while making 1~\ smaller jy,0]. The 
large experimental and theoretical uncertainties in Tj make testing of the mixing hypothesis 
via radiative transitions difficult. 

As evident from Table 13 the naive non-relativistic calculations tend to overestimate 
absolute values of the transition rates. Relativistic jf| |^ or coupled- channel corrections 
are necessary for quantitative agreement with the data. The latter is not surprising since 
non-relativistic calculations also overestimate i^{2S) — > JXcj transition rates jl6|. 
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